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Abstract 

In this paper we consider a n-dimensional stochastic differential equation driven 
by a fractional Brownian motion with Hurst parameter H > 1/3. After solving this 
equation in a rather elementary way, following the approach of [10] , we show how to 
obtain an expansion for E[f{Xt)] in terms of t, where X denotes the solution to the 
SDE and / : M" ^ M is a regular function. With respect to [2], where the same kind 
of problem is considered, we try an improvement in three different directions: we are 
able to take a drift into account in the equation, we parametrize our expansion with 
trees (which makes it easier to use), and we obtain a sharp control of the remainder. 
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1 Introduction 

In this article, we study the stochastic differential equation (SDE in short) 

X^ = a+ [ a{X^)dBs+ [ b{X^)ds, t e [0,T], (1) 
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where 5 is a d- dimensional fractional Brownian motion (fBm in short) of Hurst index 
-f/'>l/3,aeM"'isa non-random initial value and a : M" ^ iZ'^'" and 6 : ^ are 
smooth functions. 

There are essentially two ways to give a sense to equation ([T]): 

1. Pathwise (Stratonovich) setting. When if > 1/2 it is now well-known that we can 
use the Young integral for integration with respect to fBm and, with this choice, 
we have existence and uniqueness of the solution for equation ([T]) in the class of 
processes having a-Holder continuous paths with 1 — H < a < H, see e.g. plj. 
When 1/4 < H < 1/2, it is still possible to give a sense to ([1]), using the rough 
path theory, which was initiated by Lyons [HI E] and applied to the fBm case by 
Coutin and Qian [6j. In this setting, we also have existence and uniqueness in an 
appropriate class of processes. Remark moreover that, by using a generalization of 
the symmetric Russo-Vallois integral (namely the Newton-Cotes integral corrected 
by a Levy area) we can obtain existence and uniqueness for ([T]) for any H G (0, 1), 
but only in dimension n = d = 1, see [15j . 

2. Skorohod setting. Skorohod stochastic equations, i.e., the integral with respect to 
fBm in ([T]) is understood in the Skorohod sense, are much more difficult to be 
solved. Indeed, until now, essentially only equations in which the noise enters lin- 
early have been considered, see e.g., [T6]. The difficulty with equations which are 
driven non-linearly by fBm is notorious: the Picard iteration technique involves 
Malliavin derivatives in such a way that the equations for estimating these deriva- 
tives cannot be closed. 

In the current paper, we will solve ([1]) by means of a variant of the rough path theory 
introduced by Gubinelli in pTOj. It is based on an algebraic structure, which turns out to 
be useful for computational purposes, but has also its own interest, and is in fact a nice 
alternative to the now classical theory of rough paths initiated by Lyons [8], [9] . Although 
SDEs of the type ([I]) have already been studied in [10], we include in this present paper 
a detailed review of the algebraic integration tools for several reasons. First of all, we 
want to show that this theory can simplify some aspects of the analysis of fractional 
equations, and we wish to give a self-contained study of these objects to illustrate this 
point. Moreover, the analysis of stochastic partial differential equations in has lead 
to some clarifications with respect to pDj, which may be worth presenting in the simpler 
finite-dimensional context. In particular, our computations will heavily rely on an Ito- 
type formula for the so-called weakly controlled processes, which is not included in jTU], 
and which will be proved here in detail. 

As an application of this theory of integration we study the asymptotic development 
with respect to t of the quantity Ptf{a) defined by 

P,/(a) = E(/(Xn), tG[0,r],aeM", /gC°°(M";M), (2) 

where X"" is the solution of ([T]). In the case H = 1/2, the Taylor expansion of the semi- 
group Pt is well studied, see, e.g. [231 122] • Recently, Baudoin and Coutin [2] studied 
the asymptotic behaviour in the case H ^ 1/2. In this article, we extend their result in 
several ways: 
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1. In P], the authors considered the particular case 6 = 0. Consequently, their formula 
contains only powers of t of the form t"^ with n E N. Due to the drift part, we 
obtain a more complicated expression containing powers of the type with 
n, m G N. 

2. In the current article, we use rooted trees in order to obtain a nice representation 
of our formula. See also ^Sj for the case H = 1/2, and [llj for an application of the 
tree expansion to the resolution of stochastic equations. 

3. In the case where H > 1/2, we obtain a series expansion (fT5|) of the operator Pf, 
which is not only valid for small times as in [2j, but for any fixed time t > 0. 

Moreover, let us note that in [2j, the authors used the rough paths theory of Lyons 
[6|,[8|,[9] in order to give a sense to ([1]). Here, as already mentioned, we use the integration 
theory initiated by Gubinelli [TD], which allows a self-contained and hopefully a little 
simpler version of the essential results contained in the usual theory of integration of 
rough signals. 

There are several reasons which motivate the study of the family of operators {Pt, t > 
0). For instance, the knowledge of Ptf{a) for a sufficiently large class of functions / 
characterizes the law of the random variable X^. Moreover, the knowledge of Ptf{a) helps, 
e.g., also in finding good sample designs for the reconstruction of fractional diffusions, see 

The paper is organized as follows. In Section 2, we state the two main results of 
this paper. In section 3, the basic setup of [lOj with the aim of having a self-contained 
introduction to the topic is recalled. In section 4, we recall some facts on the Malliavin 
calculus for fractional Brownian motion and some properties of stochastic differential 
equations driven by a fractional Brownian motion with Hurst parameter H > 1/2. Finally, 
we give the missing proofs in section 5. 

2 Main results 

Before getting into a detailed description of the results contained in this article, let us 
first recall the main properties of a fractional Brownian motion (fBm in short). A d- 
dimensional fBm with Hurst parameter if is a centered Gaussian process, which can be 
written as 

B={Bt = {Bl...,Bty,t>0}, (3) 

where B^,...,B'^ are d independent one-dimensional fBm, i.e., each B^ is a centered 
Gaussian process with continuous sample paths and covariance function 

RH{t,s) = \{s''' + e''-\t-sr) (4) 

for ? = 1, . . . , (i. Recall also that B^ can be represented in the following way: there exists 
a standard Brownian motion such that we have 

Bl= f K{t,s)dWl 
Jo 
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for any t > 0, where K is the kernel given by 



K{t,s) 



ch 



S\l/2-H 



du 



for a constant ch which can be expressed in terms of the Gamma function. Moreover, the 
fBm verifies the following two important properties: 



(scaling) For any c > 0, -B*-'^-' = c^B./c is a fBm, 
(stationarity) For any h > 0, B.^h — Bh is a fBm. 



(5) 
(6) 



2.1 Existence and uniqueness of the solution of fractional SDEs 

As mentioned in the introduction, we will use for the integration with respect to fBm the 
integration theory developed by Gubinelli on which we try to give here a simplified 
overview. To this purpose, will denote by ^C^'" the space of linear operators from to 
M", i.e., the space of matrices of M"^'^. 

The results for fBm we will obtain in section 4 can be summarized as follows: 

Theorem 2.1. Let B be a d- dimensional fractional Brownian motion with Hurst param- 
eter H > 1/3 and a G M". Let 6 : R" — > M" and a : M" ^ C^'" be twice continuously 
differentiate and assume moreover that a and h are bounded together with their deriva- 
tives. Then the stochastic differential equation 

Xt = a+ f a{X:) dB, + /* h{X:) ds, for t e [0, T], (7) 
Jo Jo 

admits a unique solution in Q^^ai^"^) (see Definition \3.8\ below) for any n < H such that 
2k + H > I, where the integral a{X^) dBg has to be understood in the pathwise sense 
of Proposition ISTWl Moreover, if f & C^(R";]R) is bounded together with its derivatives, 
then f{X^) can be decomposed as 

f{Xn = f{a) + r Vf{X:) b{X:) ds + f V/(X;) a{Xt) dB„ (8) 
Jo Jo 

forte[0,T]. 

It is important to note that one of the main differences between our approach and the 
one developed in [6l E] is that the latter heavily relies on the almost sure approximation 
of 5 by a sequence {-B"; n > 1} of piecewise linear C^-processes, while in our setting this 
discretization procedure is only present for the construction of the so-called fundamental 
map A (see Proposition 13.21 below). 

2.2 Rooted trees and their application to the expansion of Pt 

To state the next main results we need to recall some properties of stochastic rooted trees, 
which have been introduced in [23| in the case of standard Brownian motion. 
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2.2.1 Recalls on rooted trees 



Definition 2.2. A monotonically labelled S-tree (stochastic tree) t with I = l{t) G N 
nodes is a pair of maps t= (i', i") 

f:{2,...,/}-^{l,. ..,/-!} 
f :{!,.. .,1}-^ A 

with A = {7, To, Tj^, G N} where jk is a variable index with jk G {1, . . . , d}, such that 
1f{i) < i, f {l) = 7 and f{i) E A \ {-f} for i = 2, . . . ,1. Let ITS denote the set of all 
monotonically labelled S-trees. 

We use the following notation: 

rf(t) = |{z:t"(z)=ro}| 

s{t) = \{z : t"W = r,„ J, ^ 0}| = /(t) - d{t) - 1 
p{t) = H s{t) + d{t) 

with p(7) = 0. 

In the following we denote by LTS{S) C LTS, where (S) stands for Stratonovich, the 
subset 

LTS{S) = {te LTS : s(t) = 2k, ke No} (9) 

with No = N U {0} containing all trees having an even number of stochastic nodes. 

Every monotonically labelled S-tree t can be represented as a graph, whose nodes are 
elements of {1, . . . , /(t)} and whose arcs are the pairs (t'{i),i) for i = 2, . . . , /(t). Here, 
t' defines a father son relation between the nodes, i.e., t'{i) is the father of the son 
i. Further, 7 = ® denotes the root, tq = # is a deterministic node and Tj-^, = Qj^ a 
stochastic node. Here, we have to point out that each tree t G LTS depends on the 
variable indices ji, . . . ,js{t) £ {1, • • • , dy^^\ although this is not mentioned explicitly if 
we shortly write t for the tree. 




Figure 1: Some monotonically labelled trees in LTS. 



Definition 2.3. If ti, . . . , tk are coloured trees then we denote by 

(ti,...,tk), [ti,...,tk] and {ti,...,tk}j 

the tree in which ti, . . . ,tk are each joined by a single branch to 0, # andQj, respectively 
(see also FigureW^. 
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(tl,...,tfc) [tl,...,tfc] {tl,...,tfc}j 

Figure 2: Writing a coloured S-tree with brackets. 



Therefore proceeding recursively, for the two examples t/ and t// in Figure [H we obtain 
t. = i[0%?,Oiy = i\rff,rD^ and tn = ({•^0|}|J^ = (W,r3}|J^ 

For every rooted tree t G LTS, there exists a corresponding elementary differential 
which is a direct generalization of the differential in the deterministic case, see also [23j. 
The elementary differential is defined recursively for some a; G M" by 

F(7)(x) = fix), F(ro)(x) = F(r,-)(x) = a^ix), 

for single nodes and by 

r • (F(ti)(x), . . . , F{tk){x)) for t = (ti, . . . , tfc) 
F(t)(a;) = i • (F(ti)(x), . . . , F(tfc)(x)) for t = [ti, . . . , t^] (10) 

[a^^'\x) ■ (F(ti)(x), . . . , F(tfc)(x)) for t = {ti, . . . , tj,- 

for a tree t with more than one node and with = (o"*'"')i<i<„ denoting the jth column 
of the diffusion matrix a. Here f^^\ h^^'' and cx-'^'^'' define a symmetric /c-linear differential 
operator, and one can choose the sequence of labelled S-trees ti, . . . ,tfc in an arbitrary 
order. For example, the /th component of h'^^'^ ■ (F(ti), . . . , F{tk)) can be written as 

Ji,..., jfe— 1 

where the components of vectors are denoted by superscript indices, which are chosen as 
capitals. As a result of this we get for t/ and t// the elementary differentials 

Jl,J2 = l \Ki=l 

n Pit ( " f)2 Ji ji 

F(t„) ^ /v">. ^'')) - E It E • 

Next, we assign recursively to every t G LTS a multiple stochastic integral by 

fj(A7) ift"(((t))=7 
^*'''™*'"''°\{'x.-(9(A-:)),.,.<iBi ift"(/(t))=., 

for < j < c/ with dB^ = ds. Here, t— denotes the tree which is obtained from t by 
removing the last node with label /(t). 
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2.2.2 Expansion of Pt with respect to time t 

We will denote by C^(M"';M) the space of all infinitely diff'erentiable functions g G 
C°°(M"';M), which are bounded together with their derivatives. Moreover set Am = 
{0, 1, . . . , d}"^ for m G N, and define the differential operators and T>^ as 

^'-E^'i, and V'^Y.^'-i, (12) 

k=l k=l 

for j = 1, . . . , rf. Finally, set V = V"^ . . . T'"'" for a muhi-index a e Am- 

Recall that the family of operators {Pt, t E [0,T]) has been defined by ([2]). To give 
the expansion of Pt we will use the following assumptions on the function / : R"' — »• M, 
the drift vector b = (&*)j=i,...,n and the diffusion matrix a = {(j'''^)i=i^...,n,j=i,...d- 

(A) We have /, b\ a''^ G C~(M"; R) ior i = 1, . . . ,n, j = 1, . . . , d. 

The following theorem gives the expression of the expansion of Pt with respect to t: 
Theorem 2.4. 1. If H > 1/3 and assumption (A) is satisfied, we have that 

d 

Ptfia) = Yl E ^(*)(«) E(2:t(l)o,i) t'^'^ + 0(t("'+^)^), ast^O (13) 

teLTS{S)ju-,js(t)='i- 
l(t)<m+l 

for any m G No- 

2. Let H > 1/2 and assumption (A) be satisfied. Moreover, assume that there exist 
M G N and constants K > 0, k G [0, 1/2) such that 



sup sup \V"f{x) \ + max sup sup 

a£Am a;GM" i=l,...,n a£Am x£R" 

for all m > M. Then we have 

d 



< ^"(771!)" (14) 



P*/(«)= E E F{t){a)E{Ul)o,i)t^^'\ tG[0,T]. (15) 

t&LTSiS)ji,...,js(t)=l 



Remark 2.5. (1) Here, note that each tree t G LTS{S) comprehends the variable indices 
Ji) • • • ; Js(t) which can take the values 1, . . . ,d although these variables are not mentioned 
explicitly by writing shortly t for the whole tree. The variables ji, ■ ■ ■ ,js{t) correspond to 
the components of the driving fractional Brownian motion and appear in the second sum 
in the formulas ffT3l) and ffTSl) as well as in each tree t of the summands. 



(2) In the case where H > 1/2, the boundedness of the coefficients is not needed for 
existence and uniqueness of the solution, see |T^. 
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(3) Although the additional assumption fll4p seems to be quite restrictive, it is however 
natural in a certain sense. Indeed, consider the trivial one-dimensional equation 

dX^ = dBt, tG[0,T], Xo = a 

for H > 1/2. Then we have clearly 

X^ = a + Bt, te[0,T]. 

By the first part of Theorem 12.41 and Proposition 15. 41 we have for this equation the expan- 
sion 

m 

E/(Xn = f{a) + ^Cfc/W(a)t^'= + ©(t^^'^+i)) 
fc=i 

with 

E(5i)'= r 0, if A; is odd, 

(^k = 

Hence the series 



k=l 

converges absolutely for any t G [0,T], for instance if there exists constants K > Q and 
K G [0, 1 /2) such that 

|/("^)(a)| < K'^^mlY 
for all m G N. But if, e.g., we have that 

with £ > 0, then we have 

oo 

fc=i 

for any t G (0,T]. Thus the condition (fT^ we require for the control of the remainder is 
quite natural, since the coefficients of the expansion have to satisfy a similar condition, 
as illustrated in this example. Furthermore similar growth conditions on the remainder 
or the coefficients are also usual in the case H = 1/2, i.e., for the asymptotic expansion of 
Ito stochastic differential equations respectively their functionals. Compare, e.g., [3] and 
chapter 5 in [7]. 

(4) In order to solve equation ([7]) and to bound the Malliavin derivative of the solution 
in the case H > 1/2, we need only a boundedness condition on the first two derivatives 
of b and a. To avoid too many technicalities, we have assumed in (fT^ that all derivatives 
P"/ are uniformly bounded in x G M". However, thanks to Proposition 14. II part (b), this 
condition could be relaxed, and we could allow a bound of the form 



sup + max sup 



OXi 



< K'^imlYil + \\x\ 



for a given g > 0. 
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For a better understanding of the previous results and as an example we consider 
SDE ([T]) in the case of tt, > 1, ci > 1 and give an expansion of Pt/(a) for m = 2. Here, we 
have to consider the trees with /(t) < 3 which are ti = 7^, t2 = {cTjJ^, ^3 = (t^)^, ^4 = 
= iWl}l)\ te = {[airy, tr = {{t^}1)\ ts = {r\al)\ t, = {r\al)\ 
tio = (r^, T^Y and tn = ([r^]^)^. However, only trees in LTS{S) with an even number of 
stochastic nodes have to be included since we have E(Xt(l)o,i) = for t G LTS\LTS{S). 
Then, we obtain 

d 

Pt/(a)=F(ti)(a) + F(t3)(a)E(Jt3(l)o,i)t+ Yl F{U){a)E{It,{l)o,i)t^'' 

ii.i2=i 

d 

+ J] F(t5)(a)E(Jt,(l)o,i)t'^ + F(tio)(a)E(Jt,„(l)o,i)t' 

+ F(tn)(a) E(Jt,, (1)0,1) + 0(^3^). 
Applying now (fTU|) and ffTTj) yields 

FJ{a) = f{a) + f{b){a) E( /' ds)t+ f"{^'\cT^'){a) E( /' dBH dBf) 

Jo ^.^^^.^^1 Jo Jo 

+ V /'KV))(a)E(r r dBlldBl-)t'^ + f"{b,b){a)E{r r ds,ds)e 
+ /'(6'(6))(a)E( / r rfsirfs)t' + 0(t^^). 



which finally results in 



^ r\ r ^ d n o2 r 



Ji=i i=i Ji,J2=i 



i = l Jl,J2 = l 

32 



1 (9^ f 

Jl,J2 = l 
Jl,J2 = l 



3 Some elements of algebraic integration 

As already mentioned in the introduction, we include in this present section a detailed 
review of the algebraic integration tools contained mostly in (TUl 112] • Moreover, we will 
give an Ito's type formula for so-called weakly controlled processes. 
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3.1 Increments 



The extended pathwise integration we will deal with is based on the notion of 'increments', 
together with an elementary operator 6 acting on them. The algebraic structure they 
generate is described in [TUl [T2], but here we present directly the definitions of interest 
for us, for sake of conciseness. First of all, for an arbitrary real number T > 0, a vector 
space V and an integer k > 1 we denote by CkiV) the set of functions g : [0, T]^ V such 
that gti -4k = whenever ti = t^+i for some i < k — 1. Such a function will be called a 
{k — 1) -increment, and we will set C^,{V) = Uk>iCk{V). We can now define the announced 
elementary operator 6 on Ck{V): 

k+l 

6 : Ck{V) ^ Ck^,{V), {6g),,..,,^, = J2{-l)'-'9t,-U:-t,^,, (16) 

i=l 

where ti means that this particular argument is omitted. A fundamental property of S, 
which is easily verified, is that 66 = 0, where 66 is considered as an operator from Ck{V) 
to Ck+2{V). We will denote ZCkiV) = CkiV) f] Ker(5 and BCkiV) = CkiV) f] lm6. 

Some simple examples of actions of 6, which will be the ones we will really use through- 
out the paper, are obtained by letting g E Ci and h G C2. Then, for any s,u,t G [0,T], 
we have 

i59)st = 9t- 9s, and i6h)sut = hst - hsu - Kf (17) 

Furthermore, it is easily checked that ZCk+iiV) = BCkiV) for any A; > 1. In particular, 
the following basic property holds: 

Lemma 3.1. Let k > 1 and h G ZCk+iiV). Then there exists a (non unique) f G CkiV) 
such that h = 6f . 

Proof. This elementary proof is included in [lOj, and will be omitted here. Let us just 
mention that ftj^,,,t^. = ^0*1.. .it is a possible choice. |-| 

Observe that Lemma [3.11 implies that all the elements h G C2iV) such that 6h = 
can be written as h = 6f for some (non unique) / G CiiV). Thus we get a heuristic 
interpretation of 5|c2(y)- measures how much a given 1-increment is far from being an 
exact increment of a function, i.e., a finite difference. 

Notice that our future discussions will mainly rely on /^-increments with k < 2, for 
which we will use some analytical assumptions. Namely, we measure the size of these 
increments by Holder norms defined in the following way: for / G C2iV) let 

||/IU= sup and C^iV) = {feC2iV);\\f\\,<^}. 

s,te[o,T] \t- s\^' 

Obviously, the usual Holder spaces C^iV) will be determined in the following way: for a 
continuous function g G CiiV), we simply set 

\\9h=\S9l, (18) 
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and we will say that g G C^iV) iff \\g\\^ is finite. Notice that || ■ ||^ is only a semi-norm on 
Ci(y), but we will generally work on spaces of the type 

CUV) = {9 ■■ [0,T] ^V-go = a, \\g\\^ < 00}, (19) 

for a given a E V, on which \\g\\^ thus becomes a norm. For h G C^iV) set in the same 
way 



\\h\\^^p = sup 

s,u,te[0,T] 





hgut 




\u - s\'y\t 


— u\ 


p 



(20) 



\\h\\f, = inf 1^ /i = ^/li, < Pi < yuj , 

where the last infimum is taken over all sequences {hi G C^lV)} such that h = hi and 
for all choices of the numbers pi G (0, z). Then || ■ ||^ is easily seen to be a norm on CsiV), 
and we set 

C^iV) := {heCsiVy, \\h\\,<oc}. 

Eventually, let Cl'^iV) = U^^iCl^{V), and remark that the same kind of norms can be 
considered on the spaces ZC^iV), leading to the definition of some spaces ZC^iV) and 

zcl^iv). 

With these notations in mind the following proposition is a basic result, which belongs 
to the core of our approach to pathwise integration. Its proof may be found in a simple 
form in [12]. 

Proposition 3.2 (The A-map). There exists a unique linear map A : ZC\'^{y) Cl'^iV) 
such that 

5A = Id2c]^+(^Y) and A6 = Id^^i+^^yy 

In other words, for any h G Cl'^{V) such that 6h = there exists a unique g = A{h) G 
C^'^iy) such that 6g = h. Furthermore, for any fi > 1, the map A is continuous from 
ZC^iy) to Ct^iy) and we have 

||A/.|U<^^||/.|U, heZC^.iV). (21) 

We can now give an algorithm for a canonical decomposition of a function g G C2(V), 
whose increment 5g is smooth enough: 

Corollary 3.3. Let g G C2{y) such that 6g G Ct^{V) for fi > 1. Then, for an arbitrary 
a E V, g can be decomposed in a unique way as 

g = 5f + A6g, (22) 

where f G Ci^a{V). 
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Proof. This proof is elementary. We include it here in order to see some simple manipu- 
lations of the objects we have introduced so far. 

The existence of the decomposition is due to the following fact: if Sg G Ct^{V), then 
it belongs to Dom(A). Thus, let us set h = A6g. Then S{g — h) = 0, which means that 
g — h & ZC2, and since ZC2 = BCi, there exists an element / € Ci such that g — h = 6f. 
Hence we have obtained a decomposition of the form (l22l) . 

As far as the uniqueness of the decomposition is concerned, if p satisfy ([221), then 
^/^ = and hence they differ only by a constant. Since p are both supposed to be 
elements of Ci^a{V), where a is a fixed initial condition, we obtain = p, which proves 
our claim. 

□ 

Let us mention at this point a first link between the structures we have introduced so 
far and the problem of integration of irregular functions. 

Corollary 3.4. For any 1-increment g G C2{V) such that 5g G Cl'^ , set 5f = (Id— A5)g. 
Then 

n 

where the limit is over any partition Ufs = {to = t, . . . ,tn = s} of [t, s], whose mesh tends 
to zero. Thus, the 1-increment 6 f is the indefinite integral of the 1-increment g. 

Proof. Just consider the equation g = 5f + A6g and write 

n n n 

Su = ^9uu+i = ^i^f)uu+i + ^i^S9)uu+i 

i=0 i=0 i=0 

n 

= {6f)st + J2{A6g)uu^,. 
Then observe that, due to the fact that A5g G C^'^iV), the last sum converges to zero, p. 



3.2 Computations in 

Let us specialize now to the case V = R'^foTd>l. We will also denote by C^'^ the space of 
linear operators from R*^ to M}, i.e., the space of matrices of M'^*^ and set CkC'^'^ = Ck{C'^'^). 
Then (C*,(5) can be endowed with the following product: for g G CnC'^'^ and h G Cm(M'^) 
let gh be the element of C„+m-i(I^') defined by 

{gh)ti,...,tm + n+l = 5'tl,...,t„^tn,.-,im+n-l' ^1 5 ' ' ' , ^m+n- 1 ^ [0,7"]. (23) 

In this context, we have the following useful properties. 
Proposition 3.5. The following differentiation rules hold true: 
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1. Let g G CiC^'^ and h G Ci(R'^). Then gh G Ci(M') and 

6{gh) = 6gh + g6h. (24) 

2. Let g G CiC^'^ and h G C2(M'^). Then gh G C2(M0 anc? 

6{gh) =6gh- g6h. (25) 

5. Lei ^ G Ca/:'^-' anc? /i G Ci(M'^). T/ien ^/i G C2(M0 and 

6{gh) = 6gh + g6h. (26) 

Proof. We will just prove fl24p . the other relations being just as simple. If g, h & Ci, then 

[^i9h)],t = gtht - gshs = gs {ht - hs) + {gt - gs) ht = gs {Sh)^^ + {6g)^^ht, 
which proves our claim. p-. 



The iterated integrals of smooth functions on [0, T] are obviously particular cases of 
elements of C, which will be of interest for us. Let us recall some basic rules for these 
objects: consider / G C^C'^''' and g G Cj^(M''), where denotes the set of smooth 
functions on [0,T]. Then the integral J f dg, which will be denoted by J'{f dg), can be 
considered as an element of C|°(M^). Namely, for s,t G [0, T] we set 



JMdg) = 




The multiple integrals can also be defined in the following way: given a smooth element 
h G C^C^'^ and s, t G [0, T], we set 



Jst{hdg) = 




In particular, for G C^{R'^^), p G C^C'^^^'^^ and f G C^C^^^'^-^ the double integral 
Jstifdpdf) is defined as 

Jstifdfdf) = (^j fdfdf^ = jy,^ {fdf) dfl 

Now suppose that the nth order iterated integral of f^^^df"^ ■ ■ - dp, which is denoted by 
JiP^^dp ■ --dp), has been defined for p G C^C^^-i'^j^ xhen, if p G C^iR"^^), we set 

jst{r^'dr---dpdp) = l^suir^'dr-'-dp) dp, (27) 

which recursively defines the iterated integrals of smooth functions. Observe that a nth 
order integral J{df'^ ■ ■ ■ dpdf^) could be defined along the same lines. 

The following relations between multiple integrals and the operator 6 will also be 
useful: 
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Proposition 3.6. Let f G C^C^^^ and g e C^{R'^). Then it holds that 

6g = J{dg), 6 {J{fdg)) = 0, 6 {J{dfdg)) = -{6f){6g) = -J{df)J{dg), 

and 

n-l 

5 {j{dr ■ ■ ■ df)) = -Y^j{dr... df+') J [dp ■ ■ ■ df) . 



1=1 



Proof. Here the proof is elementary again. We will just show the third of the relations. 

For s, t G [0, T] we have 

Jstidgdf) = {fu- fs)dgu = / fudgu - K,t, 

J s J s 

with Kgt = fs{gt — gs)- The first term of the right hand side is easily seen to be in ZC2- 
Thus 

5 {J{dgdf))^^, = - {5K)^^^ = - f,U - 9u], 
which gives the announced result. q 



3.3 Weakly controlled processes 

Recall that we have in mind to solve equations of the form 

dyt = a{yt)dxt, yo = a, (28) 

where t G [0,T], y is a R'-valued continuous process, a : M! —>■ C'^''- is a function, 
i.e. twice continuously differentiable and bounded together with its derivatives, x is a 
M'^-valued path and a G RMs a fixed initial condition. As usual in rough path type 
considerations, we will have to assume a priori the following hypothesis in order to handle 
equations like 



Hypothesis 3.7. The path x is -valued j-Holder with 7 > 1/3 and admits a Levy area, 
that is a process x^ = J'{dxdx) G C2^C^'^ satisfying 

(5x^ = 5x (g) 5x, I.e. [{Sx.'^)sut\ ihj) = [Sx'-]su[Sx^]ut, s,u,t e [0,T], i,j e {1, . . . ,d}. 



The solution to fizal) will then be expressed as a continuous function of the input a, a, x 
and x^. 

Let us now be more specific about the global strategy, we will use to solve equation 
fl2H|) . First of all, simple heuristic considerations show that, if the equation admits a 
solution, it should be a weakly controlled path, i.e., a process of the following form: 

Definition 3.8. Let z be a process in C]*(R^) with k < 7 and 2k + 7 > 1. We say that z 
is a weakly controlled path based on x, if Zq = a, which is a given initial condition in M.^ , 
and 6z G (M'^) can be decomposed into 

6z = C6x + r, i.e. {6z)st = Cs{Sx)st + Tst, s,te[0,T], (29) 
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with ( G CIC^'^ and r is a regular part such that r G C|'^(M'^). The space of weakly 
controlled paths will he denoted by QK,a(K'^); and a process z G QK,a{^) can he considered 
in fact as a couple {z^Q. The natural semi-norm on QK,a{^^) is given by 

withAf[g;C'^] defined by (C^j and AfiCC^iV)] = supo<,<T \Cs\v- 

Note that it is always possible to find k < 7 with 2^: + 7 > 1, since 7 > 1/3. With 
this definition at hand, we will try to solve equation (l28i) in the following way: 

1. Study the stability of QK,a(R'^) under a smooth map Lp ■.M.'' ^ W^. 

2. Define rigorously the integral / z^dxu = J{zdx) for a weakly controlled path z and 
compute its decomposition (PUI) . 

3. Solve equation fl^ in the space Qf^^ai^^) by a fixed point argument. 
In this section, we will concentrate on the first two points of this program. 

Let us first see, how smooth functions act on weakly controlled paths: 

Proposition 3.9. Let z G QK,,a{^^) with decomposition (p G C|(]R'';]R") and set 

z = (p{z), a = (p{a). Then z G QK,a(M"'), and it can be decomposed into 

6z = (6x + f, 

with 

( = Vip{z)( and r = V(p{z)r + [S{ip{z)) — Vlp{z)Sz] . 

Furthermore, 

Ar[z; Q.,a(M")] < c^,t (l + AT^f^; Q.,.(M")]) . (30) 
Proof. The algebraic part of the assertion is quite straightforward. Just write 

{6z)st = ip{zt) - ip{zs) =Vip{zs){Sz)st + 'pizt) - (p{zs) -Vip{zs){Sz)st 
= Vip{zs)Csi5x)st + Vip{zs)rst + v{zt) - ^(zs) - V(p{zs)i5z)st 

= Cs{^x)st + TsU 

which is the desired decomposition. 

In order to give an estimate for N[z] QK,a(R")], one has of course to estabhsh bounds 
for U[z]Cl{W)\, U[C] Cf£'^'"], U[C]C^C'^'''] and U[f-Cl'^{W)\. Let us focus on the last 
of these estimates, the other ones are quite similar. First notice that f = + with 

^It = ^vizsYst and rl^ = ip{zt) - ip{zs) - Vip{zs){6z)st. (31) 

Now, since Vy? is a bounded £'^'"'-valued function, we have 

X[f\ Ci^iR")] < II V<^||ooAr[r; Ci^{R% (32) 
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Moreover, 

\rl\ < ^l|VVlloo|(5;^)s*r < c^^[^; Cm'W - 

which yields 

M[f^] C2"(M")] < c^AT^fr; Cl''{^% (33) 

and thus we obtain 

M\f--Cl'm] < (1 +Ar2[r;Cr(M')]) , 
which ends the proof. |--| 

Let us now turn to the integration of weakly controlled paths, which is summarized 
in the following proposition. Notice that below we will use two additional notations. We 
will set M* for the transposition of a matrix M and denote by M ■ the inner product 
of two vectors or two matrices. 

Proposition 3.10. For a given 7 > 1/3 and k < 7, let x he a process satisfying Hypothesis 
3. 7\ Furthermore, let m G QK,b(>C^'"'^) with decomposition mo = b E C^'^ and 



{5m)st = [^^s{5x)st]* +rsu where d ClC^\ r a Cl"" . (34) 
Define z by Zq = a eR. and 

5z = m5x + n ■ x^— A(r Sx + 5ix ■ x^). (35) 

Finally, set 

J{m dx) = Sz. (36) 

Then: 

1. z is well-defined as an element of QK,a(IR)- 

2. The semi-norm of z in QK,a{^) can be estimated as 

M[z- Q.,.(M)] < c^l + T^-W[m; Q^A^^"'')]) , (37) 

for a positive constant Cx depending only on x and x^. The constant Cx can be 
hounded as follows: 

Cx < c \\x\^ + 1x^127] , for a universal constant c. 

Moreover, we have 

¥z\\^ < c,T^-W[m; Q.A'^'''% (38) 

3. It holds 

n 

Jstim dx) = lim V [nit^ {Sx)t^^u^^ + fit, ■ ^^u,u+i] (39) 
|n.,Ho .^^ 

for any < s < t < T , where the limit is taken over all partitions Ust = {s = 
to, . . . ,tn = t} of [s, t], as the mesh of the partition goes to zero. 



16 



Before going into the technical details of the proof, let us see how to recover in 
the smooth case, in order to justify our definition of the integral. (Notice however that 
( I39ll corresponds to the usual definition in the rough paths theory [8j , which gives another 
kind of justification.) 

Let us assume for the moment that a; is a smooth function and that m G C^IC^'^) admits 
the decomposition flH^ with G C'^C'^''^ and r G C^C^'^. Then J'{mdx) is well-defined, 
and we have 

f' 

ruudxu = rns[xt - Xs] + / [rUu - ms]dxu 



for s < t, respectively 

J'{m dx) = m5x + J {5m dx) . 
Let us now plug the decomposition into this expression, which yields 



J{mdx) = m5x + J{{jj,5x)* dx) + J{r dx) 

= m5x + n-x^ + J{rdx). (40) 

For the sake of clarity, let us give some details about the identity J^{{fJ,Sx)* dx) = ■ x^. 
Indeed, according to our definitions in Section 13.21 we have 



JstiilJ^5x)* dx) = / {[Sx]*^Hs) dxu = / l^s-[{Sx)su'^ dxu] = fJ^s- / [5x]s„ » = /i^ ■ x; 



for s <t, which proves the announced identity. Notice also that the terms m6x and /ix^ 
in ( l40l) are well-defined as soon as x and x^ are defined themselves. In order to push 
forward our analysis to the rough case, it remains to handle the term J'[rdx). Thanks 
to fHUj) we can write 

J7'(r dx) = Jim dx) — m5x — ^ ■ x^, 

and let us analyze this relation by applying 6 to both sides. Using the second part of 
Proposition 13.51 and the whole Proposition 13.61 yields 

6 [J{r dx)] = -6 [m 6x] - 6 [fi ■ x^] 

= —Sm Sx — Sjj, ■ + n ■ {5x ® 5x) 

= — [{fi 6x)* + r] 6x — 6fj, ■ x.^ + fj, ■ {6x ® Sx) 

= -6fi-x^ -r6x. (41) 

Assuming now that 5/i ■ x^ and r 6x are both elements of C2 with fi > 1, 6fi ■ x."^ + r 6x 
becomes an element of Dom(A), and inserting fj^Tl) into fHOj) we obtain 

6z = Jimdx) = m5x + II ■ — K{r 5x + 5/i ■ x^), 

which is the expression ( 135|1 of our Proposition 13.101 Thus fl35|) is a natural expression for 
J{mdx). 

Proof of Proposition \3.1(A We will decompose this proof in two steps. 
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Step 1: Recalling the assumption 2k + 7 > 1, let us analyze the three terms in the right 
hand side of (!35l) and show that they define an element of Q^^a such that 5z = (6x + r 
with 

( = m and f = /i ■ — A (r 5a; + 5/i ■ x^) . 

Indeed, on the one hand m G CiC^'^ and thus C = ''^ is of the desired form for an element of 
Qft.a- On the other hand, if m G 2^,6, is assumed to be bounded and since x^ G C'^CJ^''^ 
we get that /i ■ x^ G C^'{^. Along the same lines we can prove that r 5a:; G C^^'^ and 
5/i ■ x^ G C3^^'^(M). Since k + 27 > 2/t + 7 > 1, we obtain that r 5x + 5/i ■ x^ G Dom(A) 
and 

A (r5x + 5/ix2) G C2"+^(M). 

Thus we have proved that 



/X ■ x^ — A (r (5x + 5/i • x^) G C. 



and hence that z G 2^,0(1^) • The estimates (!H7j) and ([3HD are now obtained using to the 
same kind of considerations and are left to the reader for the sake of conciseness. 

Step 2: The same kind of computations as those leading to (14T|) also show that 

6 (^m 6x + fi ■ x^) = — [r 6x + 6fi ■ x^] . 
Hence equation (135!) can also be read as 

Jimdx) = [Id — A6] {mSx + ji ■ x^), 
and a direct application of Corollary 13.41 yields (!39l) , which ends our proof. 



Notice that the previous proposition has a straightforward multidimensional extension, 
which we state for further use: 



Corollary 3.11. Let x he a process satisfying Hypothesis \3 . 7| and let m G QK,fe(£'^'^) with 
decomposition uiq = b E C^'^ and 

{5m')st = [f^i{Sx),t] * + with ft' G C^C^'"^, r' G Cl^^C^^^, i = l,...,k, (42) 

where we have considered m as a 'R^^'^-valued path and have set m* = m{i, ■). Define z 
by Zq = a eM.^ and 

6z' = J{m' dx) = m'5x + -^^ + A(r^ 5x + 5^i' ■ x^). (43) 

Then the conclusions of Proposition \3.1U still hold in this context. 

Notice also that our extended pathwise integral has a nice continuity property with 
respect to the driving path x. See also [TU], p. 14]. 
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Proposition 3.12. Let x be a function satisfying Hypothesis \y. 7| and assume that there 
exists a sequence {x"; n > 1} of piecewise -functions from [0,T] to such that 

lim Ar[x" - X- C7(M'^)] = 0, and lim ATfx^'" - x^; Cl^ C^^^] = 0. (44) 

n— >oo n— +CX) 

For n > 1, define G (M'^) in the following way: set zq = h E"^^ and 

where G CIC^'^ and r" G C|''(M'^). Moreover, let z he a weakly controlled process with 
decomposition ^23\) and assume that 



lim U[z'' - z- Ci"(M'^)] + U[C - C; C^C^'^] + M[C - C; C^C^'^] + ATfr" - r; (R'^)] = 0. 

n— >oo 

Eventually, let : — a Cl -function. Then 

lim AT [jr((/?(z")rfx") - J{^{z)dx)] c:^{W)] = 0. 

n— >oo 

3.4 Stochastic calculus with respect to a rough path 

In this section, we will apply the previous considerations to two of the usual main aims in 
the theory of stochastic calculus: to study differential equations driven by a rough signal 
and to establish a change of variable formula. 

3.4.1 Rough differential equations 

Recall that we wish to solve equations of the form ( l28ll . In our algebraic setting, we will 
rephrase this as follows: we will say that y is a solution to fl28l) . if yo = a, y E QK,a(M^) 
and for any 0<s<t<Twe have 

{Sy),t = JsMy)dx), (45) 

where the integral J'{a{y) dx) has to be understood in the sense of Proposition 13. 101 Our 
existence and uniqueness result reads as follows: 



Theorem 3.13. Let x he a process satisfying Hypothesis \3. 7| and a : — C^'^ he a C 



■^2 



function, which is hounded together with its derivatives. Then 

1. Equation (^3) admits a unique solution y in QK,a(I^') for any k < 7 such that 



2ft + 7 > 1. 

2. The mapping (a, x, x^) y-^ y is continuous from R' x CjiW'-) x C^^ iW^'^'^) to 2^,0 (R^- 

Proof. We will identify the solution on a small interval [0, r] as the fixed point of the map 
r : Q«,a(ffi') ^ QkA^^) defined by r(^) = z with z = a and 5z = J{(r{z) dx). The first 
step in this direction is to show that the ball 

Bm = {z; zo = a, X[z; Q.,a(R')] < M} (46) 
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is invariant under F if r is small enough and M is large enough. However, due to Propo- 
sitions 13.91 and 13.101 and assuming r < 1 we have 

M[T{z)- Q«,J < c,,,. (1 + r^-W2[z; Q,,J) . (47) 

Since the set A = {u E : Co-,x(l + t'^^^v?) < u] is not empty as soon as r is small 
enough (see also the third point of the proof of Lemma 15731 below) . it is easily shown that 
the ball Bm defined at ( H6l) is left invariant by F for r small enough and M in A. 

Now, since we are working in Bm-, the fixed point argument for F is a standard argu- 
ment and is left to the reader. This leads to a unique solution to equation ( H5l) on a small 
interval [0,r]. One is then able to obtain the unique solution on an arbitrary interval 
[0, kr] with /c > 1 by patching solutions on [jr, (j -|- l)r]. Notice here that an important 
point, which allows us to use a constant step r, is the fact that the estimate (1471) does 
not depend on the initial condition a, due to the fact that a is bounded together with its 
derivatives. q 

Remark 3.14. The case of an equation of the form 

dyt = h{yt) dt + a{yt) dxt, yo = a, 

which can be written equivalently 

6y = J{b{y) dh) + J{cy{y) dx), yo = a, (48) 

where 6 : — is a C^-function and h : [0, T] ^ M is defined by hf = t, can be solved 
easily thanks to the previous results by taking one of the following two observations into 
account: 

1. One can define the integral J'{h{y)dh) in the usual Riemann sense for a weakly 
controlled path y. Then the fixed point argument of Theorem 13. 131 can be extended 
trivially to the case of an equation with drift. 

2. One can also define a path x = (x, /i'-'^^), where h^''^ represents k copies of the 
function h, and write equation (HSl) as 6y = J'{a{y) dx), for a new matrix- valued 
function a. The existence and uniqueness result follows then directly from Theorem 
13.131 if 6 is a C|-function. 

Remark 3.15. We have stressed here the fact that one could deal in a rather elementary 
way with processes having a regularity 7 > 1/3. However, if 7 > 1/2, our algebraic setting 
also applies and the results we have obtained so far can be expressed in a simpler way: 

Let X be a R'^- valued 7-Holder function with 7 > 1/2, and m a function in C^C^'^, 
with 7 + K > 1. Define z hj zq = a eM.'^ and 

6z = m6x — A{6m6x), (49) 

and set 

J{mdx) = 5z. (50) 

Then: 
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1. z is well-defined as an element of (K*^), and it holds: 

\J'st{indx)\ < ||m||oo||a;||^|t — + c^^K||m||K||x||^|t — sl'''^"'. (51) 

2. The integral J{mdx) coincides with the usual Young integral, and in particular, it 
holds 

n 

Jst{m dx) = lim V m*^ {Sx)ta,^^ (52) 
in., HO .^^ 

for any < s < t < T, where the limit is taken over all partitions Ugt = {s = 
to, . . . ,tn = t} of [s, t], as the mesh of the partition goes to zero. 

3. Equation (145|) can be solved in the Young sense whenever a is a C^-function, and 
the solution y satisfies 

llz/IU < Ca,a,T||a;||^, (53) 

for all A < 7. 
3.4.2 Ito's type formula 

In the sequel of this paper, it will also be essential to have a change of variable formula for 
weakly controlled process. This will be achieved under the following additional assumption 
on x^, which will be shown to be valid in the fractional Brownian motion case: 



Hypothesis 3.16. Let x be the area process defined in Hypothesis \3. 7\ and denote by 
x^'^ the symmetric part of x."^ , i.e. x^'® = ^(x^ + (x^)*). Then we assume that for 
0<s<t<Twe have 

Remark 3.17. It is worth noticing at this point that this assumption does not involve 
any limit type property of the form (jUl) for x^. This will simplify the verification of the 
different hypothesis for the fractional Brownian motion with respect to [HI [8] . 

With these assumptions in mind, our change of variable formula reads as follows: 

Proposition 3.18. Assume that x satisfies Hypothesis \3.7\ and \3.16i Let m G Qn,b{^^'^'^) 
be a process of the form (f^l) and let z G Q^^ai^^) defined by zq = a and 5z = J{mdx), 



which is given by ( f73| ). Let also f G C^(]R^;M). Then f{zo) = f{a) and 

[S{f{z))L = Jst{iV*f{z)m)dx). (54) 

Proof. The strategy of our proof is quite straightforward. By using the composition and 
integration rules for weakly controlled processes we will compute the decompositions of 
6{f{z)) and J'{(V* f{z)m) dx) respectively, and then show that they coincide. Let us 
begin with the decomposition of 6{f{z)). Recall that the decomposition fH^ of 6m can 
be written as 

Sm' = [jj,' 6xY + r\ z = 1, . . . , /c. 
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and that we have Sz'^ = J'{m^dx). Thus, for i < k the decomposition of is given by 

Sz' = m' 6x + fi' ■ - A (r* 6x + 6fJ ■ x^) . 

In the sequel of the proof, we will use the following notation: we write f for any increment 
in C2 with > 2k + 7 > 1, whose exact expression can change from line to line. In the 
same spirit, we will denote by f^'' any increment in C|'^ with regularity at least 2k. With 
these conventions in mind, some elementary algebraic manipulations yield 

1=1 *)i=i 

1=1 *ii=i 

which is the decomposition we were looking for Sf{z). 

Let us compute now the decomposition of 6[V*f{z)'m]. We have due to Proposition 
Othat 

k 

6[V*f{z)mU = Y,6[dJ{z,)l,ml + dj{zs) [6m^]^^ . (56) 

i=l 

Recall also that, setting m = [m^, . . . , m'^]*, 6z can be decomposed into 6z = m6x + f^*^. 
Thus, according to Proposition 13.91 one gets 

k 

5g{z) = ^ djg{z) Sx + f^'', 
i=i 

for any smooth function (7 : R'^ ^ M. Plugging this equality into (l56ll . we obtain that 

6[V*f{z)mU = Ast + Bst + fl^, 

where 

k k 

Ast = dy{zs)m{{5x)stml, and B^t = ^9i/(z,) [iJ^{5x)stY ■ 

i,j=l i=l 

Notice that in the definition of Ast, fn\ has been replaced by m*, since the difference 
between the two expressions is again a remainder f^'^. Now, a little elementary linear 
algebra shows that 

k 

Ast = {Sx):,Y9yizs){mirml, 
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and hence the decomposition of V*f{z)m as a weakly controlled process can be written 

as 

6[V*f{z)m] = [u5xY + r^^, (57) 

with 

k k 

= E dyiz^im^rm^ + Y,d^f{z,)fil. 

i,j=l i=l 

With the expression ( 1571) at hand, we are now ready to compute J'(V* f{z)mdx). 
Indeed, using Proposition 13.101 we get 

Jst(y*f{z)m dx) = V*f{Zs) ITLs {5x)st + l^s ■ + Tst 
k k 

= E ^^/(^^) K i^^)st + K ■ ^i) + E [i<y^i] ■ + ^st. (58) 
1=1 «j=i 



If we now put the expressions f looj) and fl58|) together, we end up with 

[6f{z)]st-Jst{V*f{z)mdx) 



1 



+ rst. (59) 



Let us show that this last expression depends only on the symmetric part x^'^ of x^. 
Indeed, it is easily checked that, if if is a symmetric matrix of M*^'*^, {M*; i <k} a. family 
of elements of M}^'^ and X e W^''^, then 



E^(^,J) [(MTM^]-X 



k d k 
i,j=l a, 13=1 * J=l 

where X® denotes the symmetric part of X. Applying this identity to H = IIess(/(xs)), 



M* = ml and X = x^^, equation (1591) becomes 



[5/(^)].t- J.i(V7(^)^^^a;) 



a 



thanks to Hypothesis 13.161 We thus have shown that 

5f{z) - J{V*f{z)mdx) = f. 



(60) 



for an increment f G with /i > 1. Now we are in the position to prove easily that 
f = 0. If we apply S to the expression above, we find that f G keiS. Thanks to Lemma 
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13. H there exists a function g & Ci such that r = Sg. Moreover, g inherits the regularity 
of f, and hence (7 G with > 1, which means that g is a constant function and that 
r = 6g = 0. Putting these considerations and equation ( I6OI) together, we finally get 

5fiz)-JiV*fiz)mdx) = 0, 

which finishes our proof. |-| 



3.5 Application to the fractional Brownian motion 

All the previous constructions rely on the specific assumptions we have made on the 
process x. In this section, we will show that the results given at Sections 13.31 and 13.41 
apply to the fractional Brownian motion. 

The combination of the following proposition with the results for the general theory 
allows us to prove Theorem 12. 1[ 

Proposition 3.19. Let B be a d- dimensional fractional Brownian motion and suppose 
H > 1/3. Then almost all sample paths of B satisfy the Hypothesis \3.7\ and \3.1(^ 

Proof. Let us first check Hypothesis 13.71 It is a classical fact that B E Cl for any 
1/3 < 7 < iJ, when i? is a fractional Brownian motion with > 1/3. As far as x^ is 
concerned, a natural choice is 

^1= dBu® dB,, i.e. x^,(z,j) = / dB^ dB^ j G {1, . . . , rf}, 

J s J s J s J s 

for < s < t < T, where the stochastic integrals are understood in the Stratonovich 
sense. Then it is a classical result that x^ is well-defined for H > 1/3 (see, e.g., [20] for 
i ^ j and [5J for i = j). The substitution formula for Stratonovich integrals also easily 
yields that 5x^ = 6x <^ 6x. Furthermore, by stationarity (EI) and the scaling property ([5]) 
of the fractional Brownian motion, we have that 

E[\^U^,J)\']=\t-s\'^E[\^l{^,J)\^]<c\t-sr. 

/^From this inequality and thanks to the fact that x^ is a process in the second chaos of 
the fractional Brownian motion B, on which all norms are equivalent for p > 1, we get 
that 

E[\^UhjW]<c,\t-s\'^^^. (61) 

In order to conclude that x^ G C'^^'C'^''^ for any 7 < 1/3, let us recall the following 
inequality from [TU]: let g G C2{V) for a given Banach space V; then, for any k, > and 
p > 1 we have 



\g\U<c{U,+2/pA9) + \\5g\U) with U,,,{g) = i / . (62) 





9st\ 


p 


\t 


— s 


7P 



By plugging inequality (I6TII into (l62ll and recalling that 5x^ = Sx ® 6x, we obtain that 
x2 G C^^C^'^ for any 7 < iJ, which shows that B satisfies Hypothesis 13.71 

The proof of Hypothesis 13.161 is now a consequence of the Ito-Stratonovich formula for 
the fractional Brownian motion (see, e.g., pL|). 
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Remark 3.20. Proposition 13.121 implies that the theory of rough paths presented here 
and the classical one of Lyons' type coincide for the fractional Brownian motion with 
> 1/3. In particular, consider the multiple integrals 

/t pui pan 

dB:i / dB:-^ ... dB:::l, for {a„ . . . , e {0, . . . , dr, 
J s J s 

with the convention B^ = t. Then these multiple integrals, which are constructed by 
means of Proposition 13.101 coincide with the usual Stratonovich integral with respect to 
the fractional Brownian motion, see, e.g., [HE]. 



4 Malliavin calculus with respect to fBm 

In this section, we assume that the Hurst index of B verifies H > 1/2. Let us give a 
few facts about the Gaussian structure of fractional Brownian motion and its Malliavin 
derivative process, following Chapter 1.2 in [TT] and Section 2 in [19j. Let £ be the set 
of step-functions on [0, T] with values in W^. Consider the Hilbert space Ti defined as the 
closure of £ with respect to the scalar product induced by 

d 

((l[o,ii], • • • , l[o,td])' (l[o,si], • • • , l[o,Sd]))^ = ^-R/f(ti, Si), Si,ti G [0,T], i = 1, . . . 

i=l 

where Rnit, s) is given by (jlj). The scalar product between two elements (f),ip & £ is given 
by 

d „i „i 

{cp,^lj)n = 7Hj2 / ^'Wiu)\r-u\''''~^drdu (63) 
Jo Jo 

with 'fH = H{2H — 1). The space 7i contains L^([0, T]; R'^) but its elements can be 
distributions, see, e.g., [21]. Formula ([63]) holds also for y?,?/; G Lh ([0, T]; M^). The 
mapping 

d 

(l[o,ii],---,l[o,td]) ^ ^Bl. 

1=1 

can be extended to an isometry between Ti and the Gaussian space Hi (B) associated with 
B = {B^, . . . , B'^). We denote this isometry by i— > B{ip). Let S be the set of smooth 
cylindrical random variables of the form 

F = f{B{^i), . . . , 5((^fc)), y., G z = 1, . . . , A;, 

where / G C°°(M'^,]R) is bounded with bounded derivatives. The derivative operator D 
of a smooth cylindrical random variable of the above form is defined as the ?i-valued 
random variable 

DF = Y,^iB{^i)^-.-^B{^^))^,. 
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This operator is closable from Lp{Q) into L^(f2;7i). As usual, D^'^ denotes the closure of 
the set of smooth random variables with respect to the norm 

||F||^2 = E|F|2 + E||L)F||^. 

In particular, if D^F designates the Malliavin derivative of a functional F G D^'^ with 
respect to B\ we have D^B^ = (5jjl[o,i] for i,j = l,...,d. 

The divergence operator 6 is the adjoint of the derivative operator. If a random 
variable u G L'^{Q]T-l) belongs to dom(5), the domain of the divergence operator, then 
S{u) is defined by the duality relationship 

E{F6{u)) = E{DF,u)n, (64) 

for every F G D^'^. Moreover, if u G dom(5) and F G D^'^ such that Fu G L'^{^; H), then 
we have the following integration by parts formula 

6{Fu) = F6{u) - {DF, u)^- (65) 

The following proposition is well known. For part (a) see, e.g., [IH] and for part (b) 
and (c), see Proposition 19 in [19] and Theorem 3.1 in [T3] . 

Proposition 4.1. Let U and a*'-', i = 1, . . . ,n, j = 1, . . . ,d be twice continuously dijfer- 
entiable with bounded derivatives. 

(a) Then equation (Q]) has a unique solution X = {X^, . . . ,X^) in the Young sense in the 
class of all processes having a-Holder continuous sample paths with 1 — H < a < H . 

(b) It holds 

max E sup \Xl\^ < oo 

i=l,...,n Q<t<T 

for all p > 1. 

(c) Moreover, we have XI G ©-"^'^(Ti) for all t G [0,T], i = l,...,n. The Malliavin 
derivative satisfies almost surely: 

DiXl = a^'^{X,) + J2 / bl^{X^)DiXtdu + J2Yl / <^u)DlXtdBl s<t, 

k=l k=l j=l 

DiXi = 0, s>t, 
for j = 1, . . . ,d, where DiXl is the j-th component of DgXl. Furthermore 

max max sup E|Z)^Xj|^ < oo. 

j=l,...,d i=l,...,n o<s<t<T 



5 Proof of Theorem 2.4 



In the present section we will prove Theorem 12.41 We separate the proof in two parts: 
in the first one, we will show how to use trees for the parametrization of the expansion; 
while in the second one we explain how to control the remainder term, which appears 
when we expand Ptf{a) with respect to t. 



26 



5.1 Rooted trees approach 



In this section, we assume that the Hurst index of the fBm B verifies H > 1/3. The first 
step in the proof of the algebraic part of Theorem 12.41 is the following result. 

Theorem 5.1. Let {X^)t(z[Q^T] be the solution of SDE ([7|j with initial value Xq = a G M". 
Then for m e Nq and f e C^+\W;R),b, , e C™+i(R";M"), 1 < j < d, we get for 
t e [0, T] the expansion 

d 

f{X1)= J2 E ^(^)(«)^^(l)o.t + 7^™(0,^) (66) 



t£LTS ji,...,j,(t)=l 
l{t)-l<m 



with a truncation term 



t&LTS ji,...,js[t)='^ 
l(t)-l=m+l 



Proof The proof is very similar to the proof of Theorem 4.2 in [23]. Recall we defined 
the differential operators and T>^ as 

k=l k=l 

for j = 1, . . . ,d. Moreover, set 

A''{[s,t]) = {(n, ...,Tk)e[o,T]'' ■.s<n<T2<---< n < t} 

for < s < t < T and A; G N. 

By reapplication of the change-of- variable formula ([8]), which holds true for the solution 
X" to our SDE, and setting = . . . V°"' and dB^'is, Si, . . . , Sk-i) = dB'^^dB'^^ . . . 
dBf^__^ for a multi-index a = (ai, . . . , ak) G {0, 1, ... , d}^, we get that 

m „ 

/(Xn = /(a) + V V Vy{a) / dB^{s,s,,...,Sk-i)+ni{0,t) (69) 

73 JAH\0,t]) 



k=l «G{0,l,---,4'' 

with the truncation term 



7^;(o,^)= Yl //•••/ '^V{x:)dBTdBT,...dB::^\ (70) 

Clearly, for m = there exists only the tree t = 7 G LTS with /(t) — 1 = and we obtain 

Y F(t)(a) Jt(l)o, = F(7)(a)=/(a). (71) 



teLT5 
Z(t)-1=0 
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Thus, to prove fl66l) it is sufficient to show for every m G N that 

^"/(«) / (i5-(s,si,...,s„_i)= Yl mio^)u^)o,t. (72) 

aG{0,l,...,4'" J^-^m]) teLTS 

l{t)-l=m 

The proof proceeds by induction. Step m = 1 is performed for a better understanding. 
In this case two different trees ti = (tq)^ and t2 = ('^fj^ in LTS, all of length 2, with 
p(ti) = 1 and p(t2) = H have to be considered. For these two trees we have 



J:e{0,l,...,d} 



k=l ji=l k=l 

F(ti)(a) Xt,(l)o,i + 5^ F(t2)(a) Xt,(l)o,t 



ji=i 



J] 5^ F(t)(a) J,(l) 



teLT5 ji,...J,(t)=l 

«(t)-l=l 



Under the assumption that equation (17^ holds for m G No we proceed to prove the case 
m + 1. Therefore, writing a = (ai, . . . , Om+i) for an element of {0, 1, ... , rf}"*"'"^, we get 



E ^"/(«) / 

EE E 'D''--'F{t){a) fui)o,s^dBZ 

&LTS ji,...,7,m=l«m+ie{0,l,...,d> 



dB°'(s, Si, . . . ,Sr 

ae{0,l,...,d}" 

d .t (73) 



teLTS ji,...,i,(t)=l am+iG{0,l,...,4 
i(t) — l=m. 

Now, we apply Lemma 2.7 in [23] to P"'"+iF(t)(a). Then, it holds for any u G LTS* with 
/(u) — 1 = m in the case of Om+i = that 

d n d d 

E ^HW = E^'(«)5^ E ^HW= E E ^WW' 

iivjs(u)=i fc=i iivJs{u)=i te-D(u) ii,.--,is{t)=i 

(74) 

where -D(u) is the set of all trees t G LTS with Z(t) = m + 2, t'|{2,...,m+i} = u', 
t"|{i,...,m+i} = u" and t"(m + 2) = tq. Clearly s(u) = s(t) holds for all t G D{vl). 

Then we proceed by considering the case of a-m+i G {1, . . . Again, by applying 
Lemma 2.7 in [23], we get for u G LTS* with Z(u) — 1 = m that 

E ^""""^ E ^(")(«)= E E^''"'""H«)^^(")(«) 

am + ie{i-,--;d} il.---j3(u)=l ilvJs(u)i"m + l=l fc = l 

= E E 

te5(u)ii,.--Js{t)=i 
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where S'(u) denotes the set of trees t G LTS with l{t) = m + 2, t'|{2,...,m+i} = u', 
t"|{i,...,m+i} = u" and t"(m + 2) = Tj_^^^^_^^. Here we have s(t) = s(u) + 1 for all t G ^(u). 

Combining now the results for the case of Om+i = and Om+i £ {^,---:d}, the 
equation 

E E ^HW= E E (76) 

«m+ie{o,i,...,4 ii,---js(u)=i tez)(u)u5(u) ii,...,is(t)=i 

holds for every u G LTS* with /(u) — 1 = m. Now it is easily seen that 



y D{u) U 5(u) = {t G LT5 : /(t) - 1 = m + 1}. (77) 

ueLT5 
/(u) — l=m 

As a last step, we observe that 



*X„(l)o,.„rfSrr^ =Jt(l)o, (78) 



holds for all t G -D(u) in the case of ctm+i = and for all t G 5'(u) with t"(m + 1) = Tq,^^-^ 
for all Om+i G {1, . . . , (i}. By applying (!7U|) - (!75|) to (I7HI) we thus arrive at (1721) with m 
replaced by m + 1, which completes the proof of the first part of Theorem 15.11 

Finally, we have to prove that 7^^(0, t) = 71*^(0, t) given in ([HSD and dnHD- From 
(jTTll it follows that for each a G {0, 1, ... , rf}™"*"^ there exists a subset LTS (a) C LTS 
with a fixed choice of ji, . . . , js(t) G {1, . . . , rf} for t G LTS{a) such that /(t) — 1 = m + 1, 
t"(l) = 7, t"(i) = r„^_, for i = 2, . . . , m + 2 and 

P"/(a)= E (79) 

tgLT^Ca) 

for all a G and m G Nq. Then, the sets LTS{a), a G {0, 1, . . . , ci}'"+\ build a 
partition of {t G LTS : /(t) — 1 = m + 1}. Thus, we have Xt(Z")o,t = Xu(^s)o,i for all 
t,u G LTS (a) and any integrable process Z. Replacing now by F(t){X^) yields that 
for all a G {0, 1, ... , d}"^~^^, the following relation holds true: 

E Xt(F(t)(x:))o,= / / .../ v'^f{x:)dB:^dB:^\..dB:-^^ (so) 

which completes the proof. □ 
Invoking Theorem 15.11 we obtain the following corollary: 

Corollary 5.2. Let {X^)tei be the solution of SDE (Qp with initial value Xq = a G M". 
Then for m e No and f e C™+^(M", R), 6, a^' G Cfe"*+^(M", M"), 1 < j < rf, t/;e ^ei for 
t G [0, T] the expansion 

d 

P*/(«)= E E ^(^)(«)E(XKl)o,l)^''(nE(7^^(o,^)). (8i) 

teLTS(5) ji,.-,isw=l 
i(t)<m+l 
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Proof. Apply Theorem 15.11 and take the expectation in formula ( !66ll . Using the notation 
of the proof of Theorem 15.11 we observe that for each t G LTS there exists an a G 
{0, 1, . . . , d}"^ with m = Z(t) — 1 such that t e LTS (a). Then, due to ( |78l) and the scaling 
property (0) it follows that 

E(Jt(l)o,t) =E(Jt(l)o,i)t^l"'+'"-l°l = E(Jt(l)o,i) (82) 

holds with \a\ = J2iLi ^{oij^o} since we have p(t) = H \ a\ + m — \a\. q 

The sequel of the paper is now devoted to derive the announced controls on the remain- 
der term E(7lm{0, t)) appearing in ( ISTl) . according to the value of H and the assumptions 
on /, b and a. These estimates will imply easily our Theorem 12.41 



5.2 Study of the remainder term for 1/3 < H < 1/2 

We assume in this section that 1/3 < H < 1/2 and that assumption (A) holds true. Then 
we will show that for fixed m G N, we have 



E(7^„(0,t)) = 0(t 



^3) 



as t — > 0, a fact which trivially yields (fT5l) in Theorem 12.41 Furthermore, notice that the 
control fl83|) is a direct consequence of the following: 

Lemma 5.3. Let g G C^(M") be bounded together with its derivatives and X be the unique 



so 



lution to ^ in QK,a(IR") with k G 



'l -H 
2 ' 



H). For any ai, . . . ,ar & {0, ■ ■ ■ ,d} we have 



where la I 



g{X)dB''^ ...dB""' 



0{i 



r-\a\{l-H) 



), as t ^ 0, 



Proof. The more difficult setting holds when aj ^ for all 1 < J < r, that is when 
r = For this reason we will only prove the assertion in this case. Moreover, we split 
the proof in several steps. 



Step 1: Scaling. For j G {1, . . . ,r} and c > set Bu^'^^^ 
the solution of ([7]), where B is replaced by B^^\ For fixed t, we have 



5"^, and let X^"'^ denote 

u/c 



L 



A-dO,*]) 



g{X)dB''- ...dB 



"1 



£ dBl^ £^ dBl^ . . . dBtg{X,^) 
J^' dBtl £^ dBtl . . . dBtlAX,,^) 
lo dBtlf C dBt^f ■ ■ ■ It' dBt-^'^gixil) 



Consequently, in order to obtain ( |8^ . it suffices to prove that 



sup E 

te[o,T] 



dBl' 



dBl' 



dBtg{xl] 



(t) 

tr 



< oo. 
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Step 2: Fix t and set 

z,= r dBt;g{X^l), [0,1]. 
Jo 

By dSZD and dSUD we have 

Here, cs > 1 is the random constant appearing in ( 1371) (see also (l38l)). whose value will 
not change from line to line, while Cg denotes a non-random constant depending only on 
(yf, whose value can change from one line to another. Set now 

qs = f dBt::^ t'^ dBt:g{xfl) = f dBt;:^ Zt^_,, s G [0, 1]. 
Jo Jo Jo 

Similarly, we have 

By induction, we easily deduce that 

M (^^ dBt: £ dBl' . . . dBt:g{xtl). Q.,o) < CB'-Cgil+U^xfl 

Since we have \zi\ < < A/'(z, Qk,o) for a path z starting from 0, we deduce from the 
Cauchy-Schwarz inequality that fl5^ is in fact a consequence of showing 



E(cb2-) < oo (85) 

and 



sup E 

ie[o,T] 



M\Xt\ = sup E \N\Xt., < oo. (86) 

te[o,T] 



Step 3: Using that B has moments of all order and flU^ . we easily obtain by 0381) that 
(!H3|) is verified. So, let us concentrate on (ISUj) . which is more difficult. We will in fact 
only prove that E|A/'^(X, Q^^a)! < +oo, since we can obtain the control (!86l) in a similar 
way. Recall from the proof of Theorem 13.13! that X defined on [0, r] belongs by definition 
to the ball Bm given by (H6l) . where M and r verify 

M>c^,b{1 + t'~^M''). 

For fixed r, the inequality u > Co-,b(1 + t'^~'^u'^) admits solutions u iff Co^b'"^ — 4:T^~'^ > 0, 
i.e., iff T^~'^ < (4 Co-.B^)""*^. In this case, the solutions are u G [M_,M+], whereby 



M± = ^- — ■ . 

By choosing for instance 



t^-'' = {8c^,b')-' (87) 
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we obtain that 

Ar(X|j„_^, , Q,,,) < (4 - 2v^)c,,B. (88) 

Furthermore, as explained at the end of the proof of Theorem 13.131 and due to the crucial 
fact that a and its derivatives are bounded, we can in fact choose the same M for the 
bound of 5X on [r, 2r], [2r, 3r], etc. Using the triangle inequality we deduce: 

< {[Tt-^\ + 1)M. 

In other words, we deduce M{X, Q^.a) < est Co-,b^^~, see (jSTl) and (j88ll . Thus it follows 
easily that E|A/'^(X, Q^.a)] < +oo and the proof of Lemma 15731 is finished. q 



5.3 Some properties of iterated integrals in the case H > 1/2 

Let us say first a few words about the strategy we have adopted in order to get equation 
( ITSll : the key point will be again to get an accurate bound for E[7lm{0,t)], and thus we 
use estimates based on Malliavin calculus tools and explicit computations of moments for 
multiple iterated integrals with respect to the fractional Brownian motion. For a proposed 
pathwise control on the remainder 7lm{0,t) , see., e.g. in [HI Remark 7.4]. 

Before we turn to the control of the remainder in the case H > 1/2, we will establish 
first some properties of iterated integrals with respect to fractional Brownian motion. To 
do this, we require some additional notations. 

For a multi-index a G {0, 1, ... , d}^ with A; G N denote by l{a) the length of a, i.e., 
l{a) = k. Moreover set = {0, 1, . . . , d}'' for A; G N, i.e., Ak is the set of all multi-indices 
of length k. Furthermore, define for a G Ak the sets 

X = {j = 1, • • • , A; : 7^ 0}, and = {j = 1, . . . , k : aj = i}, 

for i = 1, . . . ,d and |a| = \Za\- Finally for a multi-index a G Ak and j = 1, . . . , /c we 
denote 



1 



Recall that for m G N and < ti < t2 < T we set 

A™([ti,t2]) = {(ri, . . . ,r„) G [0,T]™ : ti < n < rs < • • ■ < r„ < ta}- 
Moreover, we will use the notation 

'A''([tl,t2]) Jtx Jti 



L 



dB''= I I ■ I dB^' dB"^^ . . . dB'^^ 



for a G Ak- 

With these notations in hand, the following proposition is shown easily, and its proof 
will be omitted here. Indeed, part (a) follows immediately by the symmetry of fractional 
Brownian motion and part (b) can be shown analogously to Theorem 11 in [2j. 
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Proposition 5.4. Let k eN and a E Ak- 

(a) If \a\ is odd, then we have 



E [ dB''^ 0. 



(b) If \a\ is even, then it holds 



Lnm) (l«l/2)! ^ ^ ' ^ 



with 



V(s,«)= / 

^0 



\a\/2 



n ^«,(2i_i),a^(2ol^s(20 - ts{2l-l)\'^" "^dti . . . dtk, 

0<ti<...<tk<l 

where ©3^ is the group of all permutations of the set Za, 1h — II{2H — 1) and Sij is 

Kronecker's symbol, 
(c) It holds 



E 



Afc([0,l]) 



(7///2) 



\a\ 



see' a 



with 
W(s,a) 



/ / 

J 0<ti<...<tk<l Jo<tk+i< 



\a\ 



^^5(20 — is{2i-i)\'^^ ^ dtk+i • ■ ■ dt2kdti ■ ■ ■ dtk, 



"s{2i-l)i"6{2i) 



'0<tl<...<tfc<l ^0<tfc+i<...<t2fc<l i^i 

where denotes the group of all permutations of the set 

5a = {j = 1> • • • ' 2/1; : j e 5a orj-ke 5a}. 

Notice that part (a) and (b) of the above proposition yield a representation for the 
coefficients E(Xt(l)o,i), since 



E(Xt(l)o,i) = E / 



dB"" 



A'(t)([o,i]) 



with t"{i) = ftj e {0, 1, . . . , d}, i = 1, . . . , l{t). 



For further computations, we also need the following positivity result for iterated 
integrals of the fractional Brownian motion. 



Proposition 5.5. Let G N for i = 1, . . . ,n with n e N. Moreover let a"^* e Arm ^''^^ 
< Si < ti < T for i = 1, . . . ,n. It holds 



E 



n / dB^""' 



> 0. 
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Proof. Let t[ = k2-\ A; = 0, 1, . . . , 2'. Denote by B^'^'^^ the piecewise linear interpolation 
of B^^3 ^ with step size 2~ , i.e., 

B'r- ' = B ' + 2'(i - 4)(B ' - B '), t e li. tU. 

We have 

n „ n „ 

n / dE"'"' = hm n / d^''""' 

almost surely, due to Proposition 3.12. Since nr=i /a"'»([s t ]) belongs to a finite 
Wiener chaos, we also have 



Eff / dS"""' = limEn 



dB 



according to [i]. Note that 



A'"«([s,,i,]) 



where 



for a""" ^ and = 1 for t G [0, T]. We thus have 



= ■■■ n n '^^c^ • • • dc\ ■ ■ ■ dtT" . . . dc:. 

^A'"n([s„,t„]) J^'^^iisiM) i = l j = l ^ 

But the term nr=i Ylj^=i ^ ^ product, which consists only of increments of the 

independent fractional Brownian motions B^,...,B'^ with Hurst parameter H > 1/2 
and of the constant factors 1. Since it is well-known that the increments of a fractional 
Brownian motion of Hurst index H > 1/2 are positively correlated, and also that we have, 
for a centered Gaussian vector (Gi, . . . , G2k)- 



1 ^ 

E(G'i . . . G2k) = -jj^ W ^{G s(2t)G s{2t-i 



we clearly deduce that 



n rrii 



i=l j=l 
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for all f^\...,t2ie [0, T]. Hence we obtain 

n „ 

eTT / rfs^'""" >o 

for every / G N, and the assertion follows. 

□ 

Our estimate of the remainder will also require the Malliavin derivative of an iterated 
integral. Recall then that for a random variable F G D^'^ we denote by D^F the z-th 
component of the Malliavin derivative, i.e., DF = {D^F, . . . , D'^F). Recall moreover that 
for a G ^A;, A; G N, we have defined 

= {j = I, ■ ■ ■ ,k : aj 0}, Za,i = {j = I, . . . ,k : aj = i}, 

ioY i = 1, . . . ,d and 

a_j = (ai, ^2, . . . , Oj-i, Oj+i, . . . , Ofc). 
Then the stochastic derivative of a multiple integral can be computed as follows: 

Proposition 5.6. Let m G N and a G Am- We have 



Dl / (89) 

'A™([s,t]) 



= J2 I dB^-^{U,...,t 

fori = l,...,d. 



m— 1 J 



•c", ^ S<tl<...<t,_l <tJ<t, <...<tm-l 



Proof. We proceed by induction over /(a). 

(a) Assume that l{a) = 1. For a = (0) the assertion clearly holds. Moreover for a = (j), 
j = 1, . . . ,d, we have 

Dl f dB^J^ = DUbI^^ - = for z = l,...,d, 

J s 

which corresponds to expression (l89ll . 



(b) Now assume that ( l89l) holds for all multi-indices of length m and alH = 1, . . . , d. For 
a G ^m+i we have 



[ dB'' = D\ [ Y^dB^^^+'^ 

JA'"+i([s,t]) Js 



with 



'a^cIs.t]) 
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where a = «_(„,+!)• If «m+i = 0, then 



Dl I Y^dT= / DlY^dT. 



Hence we obtain by the induction assumption 
Di [ dB"{tu...,tm+i) 

JA'"+i([s,t]) 

= J2 f [ dB^-^{ti,...,tm-i)dT 



^ti<...<tj-i<U<tj<...<tm-i<T 



J2 / rf5"-(ti,...,t„_i,r) 



.^^ ^ S<tl<...<t,-l<U<t,<...<T<t 



which shows the assertion in this case. Hence it remains to consider the case Om+i 7^ 0. 
In this case, some standard arguments based on the linear interpolation of Y and Lemma 
1.2.3 in [nj yield 



D\ fYrdBi"-^^^= f DlYrdBi^-^^^ + YJ,,^^^,l[s,t]l 

J s J s 

for 2 = 1, . . . , (i. But now we obtain by the induction assumption that 



u) 



J2 / 



■iCZ-i- -J S<tl<...<tj^l<U<tj<...<tm~l<T 



JA™([s,u]) 



/ dB"-^{tu...,tm-l,T) 



which is our announced relation. 



/ dB"'-^ {ti, . . . ,tm-l,tm), 



□ 



Now, we will establish an estimate for the second moment of an iterated integral, 
which will be the key for the control of the remainder TZm{0, t) in the expansion of Ptf{a). 
Indeed, the term (m!)~^/^ appearing in ( l90l) will be crucial in order to get some series 
convergence which will entail a nice bound on 7^m(0,t). 
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Proposition 5.7. Let m G M and a G Am- There exists a constant Ki > 0, depending 
only on H and T, such that 



E 



1/2 



< 



f _ g\\a\H+m-\a\ 



(90) 



' A™([s,t]) 

/or < s < t < T. 

Proof. The proof is separated in tliree steps. 

(i) By stationarity ([6]) of tlie fractional Brownian motion, it follows 

c 



I dB" = [ 



dB'^. 



A™([0,t-s]) 



Hence we obtain by the scaling property ([5]) of fractional Brownian motion that 



Now from Proposition 15.41 (c) it is obvious that we have 







2 




E 




< E 






iA™([0,l]) 




iA™([0,l]) 



(91) 



where a is given by a = (ai, . . . , am) with aj = if j G Za,o and aj = 1 if j G Za, i-e., all 
integrals with respect to B^'^\ i = 2, . . . ,n are replaced by integrals with respect to B^^\ 

(ii) In the next step, we will replace also the integrals with respect to t by integrals with 
respect to B^^\ More precisely, we will show that 



E 



dB' 



A^dO,!]) 



A '"([0,1]) 



(92) 



with = H{2H — 1). To prove (|92l) assume first that there is only one integral with 
respect to t, i.e. \Za,o\ = 1- Thus we have 



dB' 



A™([0,1]) 



A''! {[0,1]) Jo J A''-2([o,s]) 



dB"' ds dB""' 



with ki + k2 + 1 = m and a = {a2, 0, ai). By rearranging the order of integration, which 
is possible, since all integrals are pathwise defined, we get 



JA'=i([0,1]) Jo JA'=2([0,s]) 



dB''^ ds dB""^ 



Ysds, 



where we have set 



dB"' dB'^K 



A'=i([s,l]) JA'=2([0,s]) 
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With this notation in hand, observe that we also have 



/ 



^5(1,...,!) = / V 



•([0,1]) 



Hence, when \Zao\ = 1, one can recast (lU^ into 



E 



Y.ds 



YsdBf^ 



(93) 



We will now proceed to the estimation of the two terms in fl95l) : first of all, we easily get 



E 



Y,ds 



Jo Jo 



Let us compute now E| YsdB'-^\s)\'^: by the relation between the Young and the diver- 
gence integral for fractional Brownian motion, see, e.g. [T] or Proposition 5.2.3 in [T7], we 
have ^ ^ ^ 

[ YJB^'\s) = 6^'\Yl[o,i])+lH [ [ DlYs,\s,-S2\^''-^ds,ds2, 
Jo Jo Jo 



where we use the notation 



<5«(FV])=<5((Fl[o,i],...,0)). 



Thus we obtain 



E 



YdB^'Hs] 



r t DlY,,\si-S2r~'dsids2 
Jo Jo 



+27^,E5«(Fl[o,i]) / / DlY,,\si-S2r-'ds,ds2 

Jo Jo 

>E|<5«(Fl[o,i])f 

+27HE5«(yi[o,i]) / / DlYs,\si-S2\^''-^dsids2. 

Jo Jo 



Since clearly D^^K-Jsi — S2p^ ^ dsi ds2 G we have, owing to flM|l . that 



E 



<5«(Fl[o,i]) / / DlY,,\s^-S2r-'dsrds2 
Jo Jo 

= 1H f f f [\[Y,,Dl^DlY,,]\s,-s^\^''-'\s,~S2\'''-'dssds,dsids2. 
Jo Jo Jo Jo 

By the definition of Ys,s G [0,1], and applying Proposition 15.61 we can decompose the 
product Ysj^Dl,^Dl,^Ys^ into a sum of products of iterated integrals, and hence 

E [Ys.Dl^Dljs,] > 0, su S2, S3, e [0, 1], 
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by Proposition 15.51 Consequently we obtain 



E 



> E|5W(ri 



[0,1]. 



Furthermore, invoking [Ij, we get 

E|5«(ri[o,i])f = 7h /' ['EY,,Y,,\si-S2r^'ds,ds2 

Jo Jo 

+ f f f f EDlY,,DU,,\s^-S2\'''-'\r,-T2\'''-'ds,ds2dT^dT2. 
Jo Jo Jo Jo 



Besides, according to Proposition I5.5[ and thanks to the fact that both Kj^Kjj 
D].^Ys^Dl^Yr2 are products of iterated integrals, we obtain that 

EYsJs, > 0, and ED\Ys,D]jr, > 0, si, S2 E [0, 1]. 



Since 
we end up with 



E / YsdB'^^\s) >7hE / Ysds 
Jo Jo 

which is the announced relation (l93l) . We have thus proved that 

2 



E 



A'=i([0,l]) ^0 JA''2{[0,s]) 



dB"^ ds dB'^' 



A'=i([0,l]) Jo JA'=2([o,s]) 



dB"' dB^^^ dB' 



Applying this procedure m — |q;| times to replace all integrals with respect to t, equation 
is now easily checked. 



(iii) Let us conclude our proof: combining (l9Tll and ( l92l) yields 



E 





2 lal 

<2^E 


L 


iA'"([0,l]) 







^5(1,...,!) 



'A™-([0,1]) 



But clearly 



and thus we have 



Since 



7a™([o,i]) 



E 



A™([0,1]) 



(2m)! 
2'"(m!)3' 



(2m)! 
2"^(m!)2 
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< 2' 



the assertion flUUl) follows. □ 



Putting together Propositions 15.61 and 15. 7[ we also get the following estimate for the 
second moment of the Malliavin derivative of an iterated integral. 



Proposition 5.8. Let m G N and a G Am- There exists a constant K2 > 0, depending 
only on H and T , such that we have 



E 



Dt I dB" 

A'"([s,t]) 



1/2 



— \\ia,i\ 



V(^-l)! 



\t-s 



{\a\-l)H+m~\a\ 



(94) 



for i = 1, . . . ,d and all < s < t < T . 

Proof. Thanks to Proposition 15.61 we have that 



^A™([s,t]) ^.g-j J S<tl<...<tj^x<U<tj<...<tm-l<t 



Thus it follows 



E 



Dl. I c/S°(ti,...,t, 

A™([s,t]) 



1/2 



3 £-»Ja ,i 



dB""-^ {ti,...,t 



m—l. 



S<ti<...<tj — i<U<tj<...<tm—l1^t 



1/2 



(95) 



Furthermore, it is easily checked that 



dB""-^ {ti,...,t 



m—l. 



S<tl<...<tj — l<U<tj<...<tm — l<:t 



dB"'' X 



A'(<»-'i)([s,u]) 



dB 



Q,J2 



with a = {a^'^,i,a^^). Since an iterated integral belongs to a finite chaos with respect to 
B, all its L'P norms are equivalent. See, e.g.. Theorem 1.4.1 in [l7]. Thus, we obtain from 
Proposition 15.71 and Holder's inequality that 



E 



dB°-^{ti,...,t 



m—l I 



S<tl<...<tj^l<U<tj<...<tm-l^t 



1/2 



C2,4 



\f — g\\a-j\H+m-l-\a. 



m—l 



^/l(^.^/l{^.' 



(96) 
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with a constant 02,4 > 0. Moreover, it is readily seen that 



< 



1 



and according to the fact that {2k)\/{k\Y < 2'^^, we end up with 



< 



2(m-l)/2 



Plugging this inequality into and we obtain 



E 



d: 



m—l, 



1/2 



(v2£l)^ 1^ I n _ |(|a_,.|//+m-l-|a_,|) 

^(m - 1)! 

and since 

\a^j\H + m — 1 — \a^j\ = {\a\ — + m — |a|, 
our claim flMl) follows. 



□ 



5.4 Study of the remainder term for H > 1/2 

To avoid notational confusion we will write in the following Xf, t G [0,T], instead of Xf, 
t e [0, T], for the solution of the SDE with Xq = a. Moreover, recall that XI, t G [0, T], 
denotes the i-th component of X. Recall also that the differential operators T>^ and T>^ 
are defined as 



P0 = ^6^ 



d 



and 



D3 



E 



0" 



d_ 



(97) 



k=l k=l 

for j = 1, . . . ,d and that we have set V"' = V"'^ . . . V^*" for a multi-index a G Ak- 

With the help of the auxiliary results contained in the previous section, we are now 
able to bound ETZm{0,t) in the following way when H > 1/2: 

Theorem 5.9. Let m G N, > 1/2 and assume that assumption (A) holds. Then there 
exists a constant K3 > 0, depending only on H , T , d and n, such that 

|E7^^(o,^)| < (u^+i+Um+iy)-^ 



for allt G [0,T], where 



Um= sup sup (E|D"/(Xi)|') 

a&Am 0<t<T 



2x1/2 



max sup sup E 



i=l,...,n 



aeAm 0<t<T 



_d_ 

dxi 



V^f{Xt 



1/2 
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and 



y = max max sup ( E I -D^X 

i=l,...,n 7=l,.--,rfn<i/<.<;<T V 



i=l,...,n j=l,...,(io<u<s<T 

Notice then that the second part of Theorem 12.41 is an immediate consequence of the 
above estimate. 

Before we can prove Theorem 15. 9[ we will need the following proposition, which is a 
straightforward consequence of Proposition 5.2.3 in [T7], Proposition 14. H and the properties 
of iterated integrals of fractional Brownian motion. 

Proposition 5.10. Let m e N, a G Am, g e Cf(]R";R), and set Ja{s,t) = X\„(j^i]) dB°'. 
Then it holds 

e(^^ g{Xs)Ja{s,t)dBl 

/ / Y1^9.XXs)Ja{s,t)DiXl\s-u\^"-^duds 

+7hE (^^ ^ g{Xs)DiJ^{s,t)\s-u\^''~^dud^ 
for t G [0, T] and j = 1, . . . ,d. 

We are now ready to prove the main result of this section. 
Proof of Theorem \5.9\ : Note that by the proof of Theorem 15.11 we have 

7^^(0,t)= J2 •••/ 1^yiX,,)dB^^^dB^^\..dB:-^\\ 



m + l 



(a) We first consider a single integrand. By interchanging the order of integration, which 
is possible since all integrals are pathwise defined, we have 



/ / ■ ■ ■ / P"/(X,J dBl^ dBl' ■ ■ ■ dBtX: 
Jo Jo Jo 

= ■■■ dBl' . . . dBt- dBtX: V'^fiXt, ) dBt,^ 

Jo Jti Jti Jti 

= [ [ dB''-'V''f{X,)dB^K 



Recall that 



and 



'0 JA"^([s,t]) 

Um= sup sup (E|P"/(Xi)|2)'/' 



Um = sup sup sup I E 

i=l,...,na&Am 0<t<T 
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2\ 1/2 



If «! = 0, we clearly have 



E / P"/(XiJrfS"(ti,...,Wi) 

A'"+i([0,t]) 



m+1 



E 



A'"([s,t]) 



1/2 



ds 



If «! 7^ we have, according to Proposition IS.lOt 



A'"+i([0,t]) 

t rs " 



(9^ 



< 



I fY. E Ai)-/(X,) I dB--^ DZ^X^ 



+1H 



'o Js JA^{[s,t]) 



Thus it follows 



JA'"+i([0,t]) 

< x^z^^+i7H r r (e|d:^x^iy^' (e / 

til Jo Jo ^ ^ \ JA'^{[s,t]) 



1/4 



Is — uP''^ duds 



So recalling that we have set 



u 



dB"- 



A'"([s,t]) 



1/2 



s — ^ (i-u (is. 



= max max sup (E|D^X*M 



we get 



A'"+i([0,t]) 



< nUm+iyiH / / I E 



+Wm+i7// / /IE 

Jo Js 



dB 



1/4 



dB"- 



s — M I duds 

Is — ds. 
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Furthermore, invoking again the equivalence of U norms for iterated integral and Propo- 
sition [5l71 we obtain 



dB 



Jo \ JA"^{[s,t]) 



1/4 



t PS 



\s — ul"^^ duds 



\t- sl^^'ls-ul^^-^duds 



m\ Jo Jo 



By Proposition 15.81 we get similarly 



7^ f f ( E DZ^ [ dB'^-- 

Jo Js \ ^A'"([s,t]) 

V{m - 1)! 

Thus, we have shown for ai ^ the estimate 



1/2 



\s — ul"^^ duds 



E / I?-/(XiJrf5-(ti,...,Wi^ 

A™+i([0,t]) 

< nWm+l3^C2,4-ft'™ 1= h Wm+iC2,4-ft'™ — , =■ 

vim - 1)! 



(b) Now we consider the complete remainder term. We have 



|E7^^(o,^)|< Yl 

aeAm+l,ai=0 



E / p"/(Xijrf5"(ti,...,Wi; 

A™+i([0,t]) 



aeAm+i,aiJ^O 



E / p"/(Xijrf5"(ti,...,Wi; 

A'"+i([0,t]) 



(99) 



Since \Am\ = {d + 1)™, it follows by flUHl) and fl^ . that there exists a constant ii's > 
depending only on H, T, n and d such that 



E7^^(0, t)| < (Um+i + 3^W„+i iT; 



^H{m+l) 



which completes the proof. 



□ 
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